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> ■ ABSTRACT 
t> 

A variety of models have been explored in regard to the origin of the excess 
extreme ultraviolet (~ 0.1 keV) emission in the Coma cluster. It is now estab- 
lished that the flux is non-thermal and the only non-thermal source mechanism 
that appears viable is inverse Compton emission produced by ~ 100 MeV elec- 
trons interacting with the cosmic microwave background photons. All but one 
of the models that have been proposed require a cluster magnetic field < lfiG. 
However, recent observations strongly suggest the magnetic field in the Coma 
cluster is ~ 5fiG. We investigate the constraints on models imposed by a 



cluster field and find a limited class of models that are compatible with this 
constraint. We also investigate the possibility that the excess hard (40-60 keV) 
X-ray emission in the cluster is produced by inverse Compton emission with the 
same electron population that produces the EUV excess. We find no scenarios 
that are compatible with a large cluster magnetic field, and consequently in this 
case these two components must be unrelated. 

Subject headings: galaxies: clusters: individual (Coma) — magnetic fields — radiation 
mechanisms: non-thermal 



1. INTRODUCTION 



EUV emission in the Coma cluster in excess of that produced by the thermal X-ray 
gas was first reported by Lieu et al. (1996). Further analysis of data on this cluster was 
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carried out by Bowyer and Berghofer (1998) and Bowyer et al. (1999). Initially, this flux was 
attributed to thermal emission from "warm" gas at 10 6 K. Maintenance of a warm intracluster 
gas is extraordinarily difficult, and on these grounds alone it was generally agreed that a 
thermal source was untenable. Observations relevant to this issue were obtained with the 
Hopkins Ultraviolet telescope (Dixon et al. 1996), and FUSE (Dixon et al. 2001a,b). No 
Far UV line emission from gas at 10 6 K was detected. More recently, observations of several 
clusters with XMM (Tamura et al. 2001; Kaastra et al. 2001; Peterson et al. 2001) have 
shown no evidence for a 10 6 temperature gas. All additional clusters examined with XMM 
also show no evidence for a "warm" 10 6 gas (S. Kahn, private communication). The sum of 
all these findings seems compelling: a thermal mechanism for the EUV excess can be ruled 
out. (However, for an alternate point of view see Mittaz et al. 1998; Lieu, Bonamente, & 
Mittaz 1999; Lieu et al. 1999; Lieu, Bonamente, & Mittaz 2000; and Bonamente, Lieu, & 
Mittaz 2001a,b.) 

Since the source mechanism is not thermal, it must be the product of a non-thermal 
process. Inverse Compton scattering (ICS) of cosmic rays with the 2.7 K cosmic microwave 
background was suggested early-on as a possible source mechanism (Hwang 1997; En/Slin 
& Biermann 1998). This is the only non-thermal mechanism that has been suggested as a 
possibility for the source of this flux. 

A number of researchers have explored models of ICS in an attempt to explain the EUV 
excess in the Coma cluster (Hwang 1997, Bowyer & Berghofer 1998, En/51in & Biermann 
1998, Sarazin 1999, Atoyan & Volk 2000, Brunetti et al. 2001). With one exception (Atoyan 
& Volk 2000), all of these models require a low (< 1/zG) field and typically these models 
require a very low magnetic field (< 0.5/iG) field. 

A hard (25-80 keV) X-ray flux in excess of that produced by the 10 8 K thermal X-ray gas 
has been detected with observations with BeppoSAX (Fusco-Femiano et al. 1999). Brunetti 
et al. (2001) have attempted to explain both the EUV and hard X-ray excesses in Coma by 
a single complex model. This model also employed a low magnetic field. 

New studies of the magnetic field in clusters of galaxies strongly suggest that magnetic 
fields in clusters of galaxies are ~ 5fiG. In this paper we investigate constraints imposed by 
a high magnetic field on models capable of producing the EUV excess in clusters of galaxies 
via the ICS process. We also briefly consider the effects of a high magnetic field on attempts 
to produce both the EUV and the hard X-ray emission by the ICS mechanism with the same 
population of cosmic ray electrons. 
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2. 



ASSUMPTIONS AND OVERALL APPROACH 



We make the following assumptions. 

A. ) The magnetic field is 5/zG. The magnetic field in the Coma cluster has been the 
subject of considerable study and observational results vary considerably. Estimates based 
in part on equipartition arguments typically result in low field estimates (Rephaeli 1988, 
Giovannini et al. 1993). However, recent observational results strongly indicate that cluster 
magnetic fields are quite large. Clarke et al. (2001) studied sixteen clusters; a substantial 
number of data points were obtained for each cluster. They found the cluster magnetic fields 
were all in the range ~ 4 — 7/zG. These results seem compelling, but it could be argued 
that they are not applicable to the Coma cluster since that cluster was not included in their 
study. However, the Coma cluster field can be determined by a different procedure albeit 
for only one location in the cluster. Feretti et al. (1995) studied the rotation measure of 
the cluster-embedded head-tail radio source NGC 4869 (5C4.81). They found the number of 
magnetic field reversals through the cluster > 200 indicating B > 4.9/xG. This measurement 
is limited by the resolution of the VLA configuration employed and, in principle, is a lower 
limit for the number of magnetic field reversals through the cluster. It is possible that with 
even higher resolution the number of magnetic field reversals through the cluster could be 
larger, with a resulting increase in the estimate of B. However, this estimate is within the 4 
to 7 fxG range found by Clarke et al. (2001) for all clusters in their sample. In view of these 
observations, we find the assumption of 5/zG for the field in the Coma cluster to be quite 
reasonable. 

B. ) A scattering process is operative in the cluster which will result in a sufficiently long 
path length for the electrons that the deposition of the energy of the electrons occurs within 
the cluster itself. We note there is overwhelming observational support for this assumption 
since if the cosmic rays did not deposit their energy within the cluster, they would escape the 
cluster and produce excess EUV emission over a region which is far larger than the cluster 



C.) We assume the evolution of the electron population is given by the equation (Sarazin 



where N^d'y is the total number of electrons in the range 7 to 7 + d'-f, 6(7) is the energy 
loss rate of an electron with energy of 7, Q(^)dj is the injection rate of electrons in the 
energy range 7 to 7 + d^. 

D.) We assume the injected cosmic-ray electrons follow a power law distribution. 



itself. 



1999): 



dNjj) _ <9[6( 7 )iV(7)] 
dt 97 



+ Q(7) 



(1) 
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E.) We assume the only energy loss mechanisms are inverse Compton emission, syn- 
chrotron radiation, Coulomb ionization and bremsstrahlung. 

The energy loss rate for relativistic electrons interacting with the cosmic microwave 
background is given by Longair (1994): 

M7) = \—1 2 Ucmb = 1.37 x lO-Vs- 1 (2) 
6 m e c 

where ot is the Thomson cross section and Ucmb is the energy density of the cosmic mi- 
crowave background. The loss rate by synchrotron emission is given by Longair (1994): 

bsynil) = \—l 2 U B = 1-30 X 10- 2 V (T^)V 1 (3) 

3m e c 1/xG 

where B is the magnetic field and Ub — B 2 /8ir is the energy density of the magnetic field. 
The loss rate for Coulomb ionization is given by Sarazin (1999): 

b Coul {l) ~ 1-2 x 10" 12 n e [1.0 + ^{^] s- 1 (4) 

id 

where n e is the thermal electron density in the ICM of the Coma cluster. The loss rate for 
the Bremsstrahlung is approximately (Sarazin 1999): 

hrem{l) ~ 1-51 X IO^VtMt) + O^S" 1 (5) 



A comparison of the energy loss rates for these four mechanisms is given in Fig.l. 
This shows that ICS is the dominant energy loss rate mechanism for relativistic electrons 
in a 0.2/zG cluster field while synchrotron radiation is the main energy loss rate process 
for relativistic electrons in 5/xG cluster field. A diagram of the energy loss rate timescale 
l/ipic + b syn + bcoui + hrem) versus electron energies is shown in Fig.2. Note that there is 
a decrease of 10% in the energy of the peak energy loss, and, more importantly, there is a 
decrease in the peak energy loss rate of ~ 40% for the case of a 5/iG field as compared to a 
0.2/iG field. 

The observed EUV ICS flux is given by Blumenthal k Gould (1970): 

Me) = ^-^ XKeX {kT) ^ X F(P) X e "^' (6) 

where the symbols have their standard meaning. The energy of the photon observed, e, has 
been obtained by folding the effective area of the EUVE Deep Survey Telescope (Sirk et al. 
1997) with absorption of EUV emission by the Galactic interstellar medium in the direction 
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of Coma following the prescription of Bowyer et al. (1999). The peak of the resultant 
telescope effective area is at 80 A. 

The synchrotron radio flux is given by Blumenthal and Gould (1970): 

i- i \ 1 4-7re 3 T ^ ^p+i , 3e . P -i . . - P +i . . 

We have used an integrated diffuse flux density for Coma C of 640 mJy at 1.4 GHz (Deiss et 
al. 1997) in combination with a spectral index a = (p — l)/2 = 1.16 (Bowyer & Berghofer 
1998) as our observational base. 

We investigate three scenarios. In the first, we consider a large injection event occurring 
over a relatively short time which produces the cosmic-ray electron population. In the second 
scenario, we assume a continuous sequence of small injection events. Finally, we consider a 
combination of these two: a large injection event followed by continuous smaller events. In 
Table 1 we show the input parameters for these three models. 

In this Table, N Q is the normalization factor and p is the spectral index of the initial 
injected rate. Q Q refers to the normalization factor and p' is the spectral index of the 
continuous injection rate. Q a is related to the injected electron number Q by the relation 
Q = Qo x 7~ p - n e is the thermal electron density in the ICM. EI denotes the energy input. 



3. RESULTS 

We first consider the case of a single primary injection event. The results with input 
parameters from Model 1 are shown in Fig.3. This population will produce the EUV emission 
as reported by Bowyer et al. (1999) shown as a cross. This population greatly overproduces 
the observed synchrotron radio emission which is shown as a heavy solid line in the figure. 
With the passage of time, synchrotron losses degrade the higher energy electrons, and the 
resultant population after 1.4 x 10 s yr and after 2.5 x 10 8 yr is shown by line A and by 
line B respectively. The results shown in this figure demonstrate why it is not possible to 
produce both the EUV and the radio emission by a population of electrons produced in a 
single initial event. 

We next consider the case of a continuous distribution of smaller events. The results 
with the input parameters from Model 2 are shown in Fig. 4. The dashed line shows the 
total injected electrons. The thin solid line shows the electron population after 2 x 10 9 yr. 
The heavy solid line shows the distribution of the electrons that are required to produce the 
observed synchrotron radio emission and the cross shows the distribution of the electrons 
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producing the observed EUV emission. The results shown in this figure show why no uniform 
continuous distribution of injection events can reproduce the observational data. 

We have independently derived a variety of electron populations using a cluster field 
of 0.2/iG; these are all capable of reproducing the EUV and radio observations. This is 
consistent with results found by a number of previous studies that demonstrated that a 
continuous distribution of injection events can produce populations that will replicate both 
the EUV and the radio data in the case of low cluster fields. 

We next consider a primary event followed by continuous injections. We use input 
parameters listed in Models 3 through 11. The population obtained using the parameters in 
Model 3 after an evolution time of 1 Gyr are shown as a dashed line in Fig. 5. We also show 
the effect of changes in the size of the initial event. These results are shown as line A and 
employ inputs from Model 4. Finally, we show the effect of changes in p. Using inputs taken 
from Model 5, we obtain the population shown as line B. 

In Fig. 6 we show the effect of changes in Q and in the spectral index of the electrons in 
the input events, p', on the overall electron population. The inputs for the dashed line are 
from Model 3. Line A shows the results of decreasing the spectral index using inputs from 
Model 6. Line B shows the results of decreasing Q D using inputs of Model 7. The results 
in this figure show that allowed variations in Q Q and p' are quite small since even modest 
changes in one of these parameters alone will result in an electron population that cannot 
produce the radio data. 

We have investigated the extent to which we can mimic the best fit electron population 
shown in Model 3 by models which change p' and Q Q (the continuous injections). Inputs are 
shown in Model 8 and 9. The resultant electron distributions are shown in Fig. 7. As can be 
seen there is only a limited amount of freedom in constructing this type of model since large 
changes in these offsetting parameters which can be used to fit the EUV data will quickly 
destroy the fit to the radio data. Good fits will be obtained if p' is ~ 2.32. 

In Fig. 8 we show the effect of different thermal electron densities on the cosmic ray 
electron population. The dashed line is the product of the inputs in Model 3; the solid 
line is the result using the inputs in Model 10. The electron densities in this model are an 
order of magnitude smaller than those in Model 3. As can be seen, the cosmic ray electron 
population is not greatly sensitive to changes in the thermal electron density in the higher 
energy range. 

In Fig. 9 we show the constraints imposed on the time of the initial event in our com- 
posite model. The input parameters are from Model 3. The dashed line corresponds to an 
evolutionary time of 1 Gyr, and the solid line is the electron population after 1.4 Gyr. 
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We also investigated the effect of a 5 /xG field with regards to the hypothesis that both 
the EUV excess and the high energy X-ray flux are produced by the ICS mechanism with 
the same population of electrons. In this work we used Models 3, 8, 9, 10, and 11 which are 
capable of reproducing both the EUV and synchrotron radio flux. We investigated whether 
any of these models were also capable of producing the high energy X-ray flux by the ICS 
mechanism. The results are shown in Fig. 10. In this figure, the top line is the observed 
excess hard energy flux (Rephaeli et al. 1999). The lower lines are the fluxes produced by 
the ICS mechanism in all five models. The flux produced by the ICS mechanism in all of 
these models is more than two orders of magnitude below the observed flux. Changes in 
the input characteristics of the injection events in these different models do not affect the 
predicted flux significantly. 



4. DISCUSSION AND CONCLUSIONS 

Atoyan & Volk (2000) are the only authors who have explored the production of an 
EUV excess in a high (> 1/xG) intracluster field. They explored a model which is different in 
formal concept from those considered here, but is similar in form to our model with a large 
primary event followed by smaller injections. In specific, they considered two populations of 
cosmic rays. One is an old, highly evolved population which has undergone re-acceleration 
and compression; this population produces the EUV emission. A second higher energy 
population of cosmic rays which has been produced more recently by secondary origin is 
invoked to produce the radio emission. While this model may be appropriate, we note 
that models invoking secondary electrons for the radio emission have not been successful in 
producing the pronounced steeping of the spectral index of the radio emission with increasing 
distance from the cluster center that has been observed in well-studied radio halos such as 
Coma and Perseus. It has not been demonstrated that the Atoyan & Volk model can 
surmount this obstacle. 

In this work we have studied general classes of models to establish which of these are 
capable of producing both the EUV and the synchrotron radio emission in a 5/xG cluster 
field. We studied models with a single initial event, models with continuous injection events, 
and combinations of these two. We find that in all models with just a single initial event, 
either the EUV emission will be underproduced or the synchrotron radio emission will be 
overproduced depending on the size of the event. With the passage of time, synchrotron 
losses (which are dominant in the high field case) will remove the higher energy components 
of the initial electron population, but the distribution required to produce the synchrotron 
radio emission is not produced. This problem is basic to all models with a single electron 
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injection event. Continuous models encounter a different problem. An electron spectrum can 
be constructed which will produce the synchrotron radio emission. However, this population 
is not capable of producing the EUV flux. 

The combination of a primary event followed by smaller continuous events can produce 
the required electron distribution. However, the initial event must have occurred at a time 
no earlier than 1.4 Gyr in the past or the higher energy electrons will have been lost and 
insufficient EUV flux will be produced. The electron energy loss rate is sufficiently steep that 
even substantially increasing the number of electrons in the initial event will not substantially 
change this time constraint. On the other hand, the time of the initial event must be no 
more recent than 1 Gyr ago or the synchrotron radio emission will be overproduced. This 
requires that the injection event must have occurred in the time span 1.4 Gyr < injection < 
1 Gyr. 

Finally, we investigated the effects of a 5/iG field if both the EUV excess and the high 
energy X-ray flux are produced by ICS mechanism with the same population of electrons. 
The results shown in Figure 10 show that the high energy X-ray emission produced via the 
ICS process falls far short of the observed flux for all of these models. We conclude it is 
impossible to produce both the EUV flux and the hard X-ray emission by ICS from the same 
population of electrons with reasonable electron injection scenarios. We note that Atoyan & 
Volk (2000) have reached similar conclusions. 

The project is supported by the National Science Council of Taiwan (grant number 
NSC 90-21 12-M-008-040) and the Graduate Institute of Astronomy, NCU. S.B. thanks the 
faculty and staff of the Graduate Institute of Astronomy, NCU, and the National Center for 
Theoretical Science of Taiwan for their support and warm hospitality during an extended 
stay. 
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Table 1. 



Model N Q p El(erg) Q Q (s x ) p' EI(erg yr 1 ) n e (cm 3 ) 



Single Primary Event 
No.l 5.5 x 10 67 2.32 1.64 x 10 61 3 x 1(T 3 



Continuous Injection 
No.2 2 x 10 49 2.32 5.96 x 10 42 2 x 10 49 2.32 1.88 x 10 50 3 x 10~ 3 



Primary Event Followed by Continuous Injection 



No.3 


9 x 10 67 


2.32 


2.68 x 10 61 


2 x 10 49 


2.32 


1.88 x 10 50 


3 x 10- 


-3 


No.4 


9 x 10 68 


2.32 


2.68 x 10 62 


2 x 10 49 


2.32 


1.88 x 10 50 


3 x 10" 


-3 


No.5 


9 x 10 67 


2.5 


6.87 x 10 60 


2 x 10 49 


2.32 


1.88 x 10 50 


3 x 10" 


-3 


No.6 


9 x 10 67 


2.32 


2.68 x 10 61 


2 x 10 49 


2.2 


4.80 x 10 50 


3 x 10- 


-3 


No.7 


9 x 10 67 


2.32 


2.68 x 10 61 


2 x 10 48 


2.32 


1.88 x 10 49 


3 x 10- 


-3 


No.8 


9 x 10 67 


2.32 


2.68 x 10 61 


8 x 10 48 


2.2 


1.92 x 10 50 


3 x 10- 


-3 


No.9 


9 x 10 67 


2.32 


2.68 x 10 61 


9 x 10 49 


2.5 


2.17 x 10 50 


3 x 10- 


-3 


No.10 


9 x 10 67 


2.32 


2.68 x 10 61 


2 x 10 49 


2.32 


1.88 x 10 50 


3 x 10- 


-4 


No.ll 


3.5 x 10 68 


2.5 


2.67 x 10 61 


2 x 10 49 


2.32 


1.88 x 10 50 


3 x 10- 


-3 
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Fig. 1. — Energy loss rates for individual mechanisms including inverse Compton scattering 
(ICS), synchrotron radiation (Syn) in two magnetic fields, Coulomb ionization (Coul), and 
bremsstrahlung radiation (Brem). EER is the "EUV-emitting electron range" and SRER is 
the "synchrotron-radio-emitting electron range" . 
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Energy Loss Rate Tirnescale 




2.0 2.5 3.0 3.5 4.0 4.5 5.0 

Log 7 

Fig. 2. — The energy loss rate tirnescale for the electrons in a 5/xG and 0.2/xG field. The 
dashed line is the tirnescale for 5/xG case and the solid line is for 0.2/zG. The thick horizontal 
lines mark the EUV-emitting electron range and the synchrotron-radio-emitting electron 
range. 
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Electron Distribution 
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Fig. 3. — The electron spectrum using the inputs in Model 1. The dashed line shows the 
distribution of the initial injected electrons. Line A shows the distribution of the electrons 
after 1.4 x 10 8 yr while line B shows the distribution after 2.5 x 10 8 yr. The cross is the 
distribution of the electrons producing the observed EUV flux and the heavy solid line is the 
distribution of the electrons producing the synchrotron radio emission. 
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Fig. 4. — The electron spectrum using the inputs in Model 2 with an evolution time of 2 x 10 9 
yr. The dashed line shows the total injected electrons. 



-16- 



Electron Distribution 
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Fig. 5. — The dashed line shows the electron distribution using inputs from Model 3 with an 
evolution time of 1 Gyr. Line A with input parameters from Model 4 shows the distribution 
of electrons with N Q increased, and line B with input parameters from Model 5 shows the 
distribution of the electrons with p increased. 
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Electron Distribution 
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Fig. 6. — The effects of changes in Q Q and p' on the electron spectrum. The dashed line is 
Model 3. Line A shows the results of decreasing p' using the input parameters of Model 6. 
Line B shows the results of decreasing Q Q using the inputs in Model 7. The evolution time 
is 1 Gyr. 
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Electron Distribution 




2.0 2.5 3.0 3.5 4.0 

Log 7 



Fig. 7. — Models mimic the best fit electron distribution (Model 3). 2.2 and 2.5 are chosen 
for p' and the resultant electron distributions are shown as the solid lines. The evolution 
time is 1 Gyr. 
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Electron Distribution 
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Fig. 8. — The effect of differing thermal electron densities, n e . The dashed line is the product 
of the inputs in Model 3; the solid line is the result using the inputs from Model 8 with n e 
smaller by a factor of 10. The evolution time is 1 Gyr. 
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Electron Distribution 
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Fig. 9. — The effect of different evolution times. The input parameters are from Model 3. 
The dashed line is the distribution after 1 Gyr and the solid line is the distribution after 1.4 
Gyr. 
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Photon Flux 
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Fig. 10. — A comparison of the observed hard X-ray flux and the ICS results of the models. 
The solid line shows the excess high energy X-ray flux. The dashed line shows the flux 
produced by the ICS mechanisms using inputs from Model 3, 10, and 11. The emission 
predicted by these three models is virtually the same. The dash-dot line shows the emission 
predicted by Model 8, and the dotted line shows the emission predicted by Model 9. 



